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1. Introduction and summary 


In earlier papers surveyed recently (1) the present author tried to separate from 
the total energy content of a metal, the structural part due to thermal lattice defects 
on the assumption that the part due to lattice vibrations was approximately described 
by a Debye function. The structural energy thus obtained was used for a study of 
the relation between the solid and liquid states. 

In the present paper we shall try a more elaborate determination of the energy 
due to thermal lattice defects in solid metals starting from the general assumption 
that the vibrational energy approaches the classical value 3R7' with increasing 
temperature, avoiding the special assumptions contained in the Debye function. 
We shall also take into account the part of the energy content due to the free electrons, 
which was not done in the earlier calculations. 

The results confirm the basic assumptions and allow a determination of the vibra- 
tional zero-point energy as well as the energy due to thermal lattice defects. The 
temperature dependence of the latter shows the existence of two different kinds of 
lattice defects. One of these which asserts itself first at high temperatures near 
the melting point is identified as lattice vacancies. Its energy of formation is of the 
order of magnitude calculated by Huntington and Seitz (2), Fumi (3), Seeger and 
Bross (4) and others, and determined by resistivity measurements on quenched 
metals by Bauerle, Klabunde and Koehler (5) and others. The other kind of defects 
represents the greater part of the structural energy. The same defects also seem 
to present themsleves in a calorimetric delay, following a rise of temperature, recently 
found by Astrém (6). Their structure, which is still unindentified, will have to fulfil 

certain conditions found by the experiments. 


ee 


; 


2. Postulate for the vibrational energy 


In order to separate from the total energy content of a metal the structural part 
ue to thermal lattice defects we need theoretical expressions for the part of the 
ergy due to the thermal vibrations of the atoms. 

Our working hypothesis is that the sum of the vibrational zeropoint energy pro 
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Fig. 1. Postulated scheme for the vibrational 
energy. 


mole #, and the thermal vibrational energy H,,, approaches the classical value 3R7' 
with increasing temperature. Thus 


E,+Eyp>3RT for T>6, (1) 


where @ is a characteristic temperature of the metal. This condition is fulfilled by 
the Debye theory, in which F,,;, is fully determined by 7 and a constant § and the 
zero-point energy is 90/8, but our assumptions are more general than those of the 
Debye theory, which are now known to be too narrow. Our hypothesis is at variance 
with those earlier theories which give a surplus energy due to the anharmonicity 
of the vibrations. Those theories, however, do not give any quantitative predictions 
for the deviations from the classical atomic heat capacity 3 and may perhaps be 
looked upon as an attempt to understand experimental results for which there was 
at that time no other explanation. 
For our purposes we express the condition (1) in the form 


B,+ Ey,» =3RT + E,, (2) 


with the condition that the ‘‘residual zero-point energy” H,, decreases from the value 
E, for T =0 to zero for 7’>6. These conditions are illustrated by the diagrams of 
fig. 1. Whether the residual zero-point energy £,, can be theoretically separated as 
due to ordered vibrations or is just a formal difference is of no importance for our 
problems. 

Our working hypothesis will be well supported by the results obtained in the 
following by its application to the measured energy content at normal practically 
zero pressure. A constant higher pressure or a pressure increasing with the temperature 
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in order to keep the volume constant will surely alter the zero-point energy and the 
temperature dependence of the residual zero-point energy, but we have no reason 
to doubt that at high temperatures the total energy of vibrations approaches the 
value 3RT independent of pressure and volume. 


3. Procedure of phenomenological analysis 


For a determination of the structural energy due to thermal lattice defects, i.e. 
defects the concentrations of which are functions of the temperature, we have to 
start from the experimentally observed thermal energy per mole 


b §e 
Exps= { CdT, (3) 
0 


where C is the atomic heat capacity at zero pressure or at a pressure of one atmosphere, 
which makes no observable difference. In a metal this observed energy is assumed to be 
the sum of three parts: the energy of the free electrons H,, the thermal part Ey, 
of the energy of lattice vibrations, and the structural energy of thermal defects E,,,, 


Eobs aq Ea a Ey in as E str (4) 


The electron energy E,; is at low temperatures given by 
Ey =" T*, (5) 


where the constant y is known for many metals from experiment. Our only chance is 
to assume that this formula is approximately valid at all temperatures. 

The total vibrational energy is the sum of the thermal part H,;, and the zero- 
point energy #,. Our working hypothesis is that this sum approaches the classical 
value 3R7' with increasing temperature. 

From equations (2) to (5) we obtain the equation 


~ 


i 
E,,— E,+ Ey={CdT-3RT-<T (6) 


0 


~~ 


with known quantities on the right side and the unknown to the left. As will be 
seen in the following the variation of this sum with temperature is such that Ep, 
#, and #,, can be determined separately. 


4. Gold 


Gold is well adapted for a first demonstration of the use of an analysis by means of 
equation (6). 
The first term to the right in this equation is obtained from calorimetric measure- 
ents. For the low-temperature range 0-300°K we have used the results of Geballe 
nd Giauque (7). For the range from 300 to 600°K the increase of energy is taken 


67 


G. BORELIUS, Energy due to thermal lattice defects in solid metals 


Table 1. Aw. Energy data in cal/mole. 


° 
K Lops 3RT r the Ezy Ey Exx Eops—cal 
0 0 0 0 351 
15 0 89 0 262 
50 70 298 0 122 
100 293 596 1 46 
150 563 894 2 17 
200 849 1192 4 5 
250 1145 1490 5 0 saat 
273 1283 1627 7 0 0 
300 1445 1788 8 1 a 
400 2054 2384 14 8 dl 
500 2675 2981 22 26 =n) 
600 3309 3577 31 55 : te) 
692 3914 4125 42 92 0 i 
905 5313 5395 71 200 3 = 10) 
1073 6470 6396 100 296 22 wai 
1232 7601 7344 132 389 88 == 
1273 7910 7588 141 413 128 <a!) 


from a compilation by Stull and Sinke (8). For the range above 600°K we have used 
the measurements by Jaeger, Rosenbohm and Bottema (9) from which we have 
taken the directly determined energy differences between the high temperatures and 
room temperature. The total thermal energy values thus obtained for a number of 
temperatures are given as H,,, in table 1. It must be pointed out that the tables of 
energy values published by many authors, which are obtained by smoothing the 
observed values using power series of the temperature, are rather dangerous in our 
case because the smoothing often conceals the real dependence of temperature. 

The third column of the table contains the values of 3R7' and the fourth column 
the thermal energy of the free electrons, calculated with the value y = 1.74 x 10-4 
cal/mole deg? found by measurements on Au at low temperatures by Corak and 
Garfunkel (10). 

The sums of the columns 2, 3 and 4, i.e. the values of the right hand side of equa- 
tion (6), are plotted against the temperature in fig. 2. The curve is easily interpreted 
in terms of the left hand side of the equation. The initial decrease of the curve means 
the gradual disappearance of the residual zero-point energy H,, in accordance with 
the scheme of the lower diagram in fig. 1. Because of the low Debye characteristic 
temperature 9 = 190 of Au, #,, has practically disappeared before the rising of the 
curve due to the structural energy sets in. The horizontal part of the curve means that 
in this interval the atomic heat capacity has its classical value 32. The position of 
the horizontal line corresponds to the energy —H,. The zero-point energy comes out 
as H, = 351 cal/mole. Taking —351 as the zero point of a new scale, to the right of the 
diagram, the curves to the left and right give the residual zero-point energy H,, 
introduced in column 5 of table 1, and the structural energy in cal/mole respectively. 

The experimental points on the high temperature side still contain the original 
accidental errors from the measurements and are now graphically smoothed by a 
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Fig. 2. Determination of structural energy af Au. 
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_ Fig. 3. Separation of the energy of the two kinds 1 
of defects in Au. 
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curve. This curve has the general shape of a solubility curve. If there were just one 
kind of thermal defects we should expect the curve to fit the formula 


a 
i. een ( ss 4 (7) 


where A is a constant and H, the heat of formation of the defects. Equation (7) 

can also be expressed as 

H, loge 1 (8) 
Ri wk 


log Egtr=log A 


which means a straight line in a log H —1/7 diagram. 

In order to check this formula we have taken the values of the structural energy 
from the curve of fig. 2 with intervals of 50 degrees from 500 to 1300°K and plotted 
their logarithms against 1/7’ in fig. 3. According to the formula this plot should 
give a straight line which within the limits of a resonable uncertainty is also the 
case up to about 900°K. This straight line corresponds to an energy 


H,= 2480 exp ( = or) cal/mole (9) 


if both empirical constants are expressed in cal/mole. Values of H, obtained from 
equation (9) are introduced in column 6 of table 1. 

Above 900°K there is, however, an appreciable difference between the total struc- 
tural energy H,,, and H,. We denote this difference #,, and have 


Ey = Egy — Ey. (10) 


The plot of log #,, against 1/7’ shown in fig. 3 gives again a straight line which means 
a temperature dependence of the form of equation (7). As in this case the constant 
A, obtained by a very uncertain extrapolation to 1/7'=0, would be uncertain even 
to its order of magnitude, we prefer to give the result in the formula 


000 /1 1 
By = 10 exp ( =~ (7 aon) cal/mole (11) 


approximately valid in the range of the measurements. From this equation calculated 
values are introduced in column 7 of table 1. 

The differences between the observed thermal energy H,,, and the sum of its 
calculated parts 


Eops—cat = Bong + My — Le — Bg pote ee Ne (12) 


are given in the last column of the table. They are well within the limits of experi- 
mental error. 

Our phenomenological analysis thus shows the existence of two different species 
of thermal lattice defects D; and D,, with the thermal energies #, and E),. 

Dy, is easily identified as lattice vacancies. The empirical energy of formation 
H, = 23000 cal/mole corresponds to 1.0 eV which within the limits of error is the 
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same as the value H,= 1.02 + 0.06 eV found by Bauerle, Klabunde and Koehler 
(5) from resistivity measurements on quenched specimens, and the value 0.95 + 0.1 
reported by Bradshaw and Pearson (11). 

The defect D, which thermally dominates over D,, in the whole range of the solid 
state has the heat of formation H, = 4530 cal/mole or 0.20 eV. The structure of this 
defect is unknown and will be discussed in section 8. 


5. Silver 


The procedure for determining the energy due to the two kinds of thermal lattice 
defects which has been described in detail in the case of gold is applicable also for 
the other metals, and we shall confine ourselves to giving the results in diagrams and 
tables. 

For Ag the thermal energy in the range below 300°K is determined from the 
values of specific heat chosen in the compilation of Kelley (12) and in the range 
above 300°K it is taken from the measurements of Jaeger, Rosenbohm, and Weenstra 
(13) on their silver with a minimum content of oxygen. The data are given in 
table 2. For the calculation of the electron energy we use the value y = 1.48 x 10-4 
cal/mole deg? found by Corak and Garfunkel (10). 

The sum of the terms to the right in equation (6) is shown in fig. 4. Also in Ag, 
with § = 215, the residual zero-point energy disappears before the structural energy 
has reached a measurable size. The zero-point energy is found to be 403 cal/mole. 

The structural energy is taken from the diagram of fig. 4 with intervals of 50 


Table 2. Ag. Energy data in cal/mole. 


°K Eops 3RT z Me By Ey | Exx Eops—cal 
0 0 0 1) 403 
20 4 119 0 288 
50 52 298 0 157 
100 252 596 1 58 
150 512 894 2 19 
200 795 1192 3 3 
250 1091 1490 5 0 aes 
273 1230 1627 6 0 
300 1395 1788 7 0 +3 
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Fig. 4. Determination of structural energy of Ag. 


degrees from 500 to 1200°K and the logarithms are plotted against 1/7’ in fig. 5. 
The analyzis of this plot has given for the first kind of defect the energy 


5160 
E,= 4370 exp ( > a) cal/mole. (13) 


The energy due to the second kind of defect, identified with the vacancies, is about 
8.9 cal/mole at 7’ = 1000 and the slope of the line drawn in fig. 5 gives the heat of 
formation H, = 27 500 cal/mole. We have thus 


vt 1 


1 
R 7p sa00)) cal/mole. (14) 


H,=27 500 cal/mole corresponds to 1.2 eV. 
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6. Copper 


Up to 300°K the experimental energy values in table 3 are based on measurements 
by Giauque and Meads (14) and from 300° to 500°on Stull and Sinke’s compilation 
of heat-capacity data (8). Above 500° we have used the results of Jaeger, Rosen- 
bohm and Bottema (15). Unfortunately these authors have not published their 
results near to the melting point which is at 1083°C. In a foot note they say that 


Table 3. Cu. Energy data in cal/mole. 


°K Eons 3RT z sie Egy Ey Exy Eps -cal 
0 0 0 0 602 
15 0 119 0 483 
50 21 298 0 325 
100 161 596 1 166 
150 383 894 2 89 
200 645 1192 3 52 
250 923 1490 5 30 
Zio 1056 1627 6 25 0 
300 1213 1788 * 19 1 
400 1807 2384 13 6 6 
500 2415 2981 20 2 24 
586 2963 3495 28 0 42 0 
693 3650 4130 38 80 0 +4 
904 5033 5388 65 183 =—6 
1072 6186 6389 92 279 22, +6 
1233 7309 7347 122 377 71 —6 
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Fig. 6. Determination of the structural energy of Cu. 


“above 960° copper evidently shows an abnormal behaviour: its specific heat then 
increases more rapidly with the temperature. The result will be discussed in a later 
paper’. We have found no such paper. The results would probably have been of 
good use for a better determination of the vacancy energy, which in fact can be 
expected to raise the specific heat rapidly in this range of temperature. For y we 
have used the value 1.60 x 10~4 cal/mole deg? found by Corak and Garfunkel (10). 

As copper has a high Debye characteristic temperature of about 0 = 325, the 
falling curve of the residual zero-point energy and the rising curve of the structural 
energy in fig. 6 overlap slightly, which gives an extra uncertainty to the estimated 
zero-point energy, H, = 602 cal/mole, of a few units. 

The energy functions for the two kinds of defects obtained from the lines in fig. 7 are 


4950 
E,=2750 exp (- 4 cal/mole (15) 
19000 /1 1 
and Ey, =12 exp (- <4 (7, = sais) cal/mole. (16) 
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Fig. 7. Separation of the energy of the two kinds 1 
of defects in Cu. 
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The energies of formation correspond to H, = 0.22 eV for the main defects and H, = 
0.8 eV for the vacancies. The latter value agrees fairly well with theoretical predic- 
tions of H, = 1 eV by Huntington and Seitz (2), 0.9 eV by Fumi (3) and 0.8 eV by 
Seeger and Bross (4). 


7. Aluminium 


For Al we have found no directly useable calorimetric determinations of the 
energy content at high temperatures. We are restricted to the compiled and smoothed 
values given by Kelley (12) and Stull and Sinke (8) shown in table 4. The results 
give therefore only a rough orientation. For y we have used the value y = 3.484 = 10-4 
cal/mole deg? found by Kok and Keesom (16). 

Al has a very high Debye temperature of about 6 = 400 which leads to an appreci- 

able overlapping of the falling and rising curves in fig. 8. Our estimation of the zero- 
point energy gives 730 cal/ mole. The analysis of the structural energy illustrated by 
fig. 9 gives the empiric formulae 


; s 4720 
EH; = 4040 exp | = ae cal/mole (17) 
13000 /1 1 
nd Ey,=7 exp ( R (7 sy) cal/mole. (18) 


he heats of formation 4720 and 13 000 cal/mole correspond to 0.20 and 0.56 eV. 
e value H,=0.56 eV for the formation of vacancies is in fair agreement with 
lues obtained from recent resistivity measurements on quenched specimens. From 
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Table 4. Al. Energy data in cal/mole. — 


ba ee ee Seal 


r Bigg 87 1) eee Egy By Hz A) plana 
0 0 0 0 730 

20 0 119 0 611 

50 12 298 0 444 
100 117 596 2 249 
150 310 894 4 142 
200 551 1192 7 82 
250 820 1490 11 49 0 
273 950 1627 13 40 1 
300 1105 1788 16 30 1 
400 1703 2384 28 10 11 
500 2330 2981 44 0 35 0 0 
600 2987 3577 63 77 2 Se 
700 3674 4173 85 136 10 0 
800 4391 4769 112 208 32 0 
900 5137 5365 141 289 77 aa) 


the extrapolated initial increments of resistivity after quenching from different 
temperatures obtained from diagrams published by Panseri, Gatto and Federighi 
(17) we obtain H,=0.5 eV. De Sorbo and Turnbull (18) report a value in excess of 
0.6 eV. Bradshaw and Pearson (19) report the value 0.76 eV. 


8. The determination by Astrém of heat of formation of lattice defects 


We have found that the temperature variations of the heat content of metals 
reveal the existence of two different thermal defects denoted D, and Dy. Dy has 
been identified with the defects which are preserved for some time after quenching 
from high temperatures and which are theoretically interpreted as vacancies. D, is 
the thermally dominating defect in the whole range of the solid state. As pointed 
out earlier by the present author the phenomenon of melting is also dominated by 
an increase of some defects other than vacancies, which evidently are the D,-defects. 
Contrary to the Dy, the D, are not, at least not to any great extent, preserved by 
quenching. On the other hand recent measurements in this laboratory by Astrém 
(6) have given the surprising result that in very pure and coarse-grained metals 
after an increase of the temperature the development of certain defects is slow enough 
to be studied by isothermal calorimetry. A discussion of different possibilities gave 
as the most probable interpretation that these defects were a part of the unindenti- 
fied defects causing the increase of structural energy discussed in earlier papers of 
the present author. We are now able to give further evidence that the defects studied 
by faGrom are the D,-defects. 

In the log H —1/T diagrams for Ag of fig. 10 the line marked £, is the same as in 
fig. 5 and the other lines are those obtained by Astrém. Astrom had found that when 
specimens of pure Ag were heated from room temperature to various temperatures 
in the range 100 to 400°C, the heat taken up in order to obtain equilibrium at the new 
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Fig. 8. Determination of the struc- 
tural energy of Al. 
P] 


; 


temperature has an exponential tail with a time constant of a few hours. The energy 
#, contained in these tails increases with increasing temperature as 


Ve 
By= Avexp (~ 74) (19) 
here H, is a heat of formation and A, a constant increasing with increasing purity 
d grain size of the specimen. In a diagram of the type of fig. 10 this gives straight 
ines. The figure shows such lines for three specimens of 99.999 % Ag with the grain 
size increasing from A to C. At least for the specimen C the slope of the line is the 
me within the limits of error as that of the line now obtained from the energy 
mtent. This means that the energy of formation is the same and the defects most 
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1 Fig. 9. Separation of the energy of the two kinds 
of defects in Al. 
2 15 1 1000/T 05 


probably identical. The same agreement exists between the lines for H, in our figures 
3 and 9, and new unpublished results by Astrém for Au and Al. The energy observed 
for Ag in Astrém’s measurements is for specimen O with a grain size of a few mm 
about 25 % of that calculated from the energy content, and for the specimen A with a 
grain size of 0.3 mm about 6%. It is imaginable that an ideal crystal would give the 
whole calculated energy and that in the actual samples only small volume parts are 
ideal. In the other parts the development of the thermal defects may be catalysed 
by other imperfections and be too rapid to be observed at least by the used calori- 
metric method. 


9. Remarks concerning the primary defects 


The structure of the thermal defects denoted D; is unknown. The theoretical 
predictions for vacancies, interstitials or combinations of these elements do not fit the 
experimental results, from which we point out two characteristic features of the 
unidentified defects. One is that the heat of formation is rather small, very near to 
0.20 eV for the four discussed metals. The other is the fact, found by Astrém, that 
the time constant t is of the order of hours for the formation of the defects in the 
undisturbed lattice and that this time constant decreases very little with increasing 
temperature. The energy of activation as derived from log t — 1/7’ diagrams is of 
the order of only 0.05 eV. The extrapolated value of t for 7’ 9 is still about 10% 
times the period ty of the atomic vibrations. The formation of a defect at a certain 
point of the lattice does not therefore imply a very great concentration of the fluctuat-_ 
ing thermal energy to this point but a very improbable combination of the movements 
of a number of atoms. There are several imaginable defects which at first sight seem 
to comply with this condition. Theoretical evaluations of the heat of formation for 
different models, which are now going on, will probably help to make the right 
choice. 
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Fig. 10. Comparison of the energy of the 
unidentified defects in Ag with the calori- 
metric results obtained by Astrém. 
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10. The zero-point energy 


The values of the zero-point energy H,, obtained from this analysis of the energy 
content, are probably fairly accurate. In table 5 these values are compared with 
probably less accurate values obtained earlier in other ways. Two years ago the pres- 
ent writer (20) evaluated the potential half U, of the zero-point energy from plots of 


Table 5. Zero-point energy in cal/mole. 


e electrical resistivity as a function of the potential energy due to vibrations and 
ttice defects. Corresponding zero-point energy-values 2 U, are inserted in the table 
gether with the differences 2U,— H, in per cent. These differences are probably 
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mainly due to the uncertainty of the determinations of 2U,. The table also contains 
the zero-point energies evaluated from the Debye theory as 9 R6/8 using the 6-values 
from the last column of the table, which are obtained from the temperature depend- 
ence of the atomic heat. In this case the differences of about 20% from the present 
values are systematic and probably due to deviations from the distribution function 
for the frequences, proportional to 7”, contained in the Debye theory. 


Stockholm, Royal Institute of Technology. 
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